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We present a method derived from Laplace transform theory that enables the evaluation of 

fractional integrals. This method is adapted and extended in a variety of ways to demonstrate its 

utility in deriving alternative representations for other classes of integrals. We also use the method 

in conjunction with several different techniques to derive many results that have not appeared in 

^ , tables of integrals. 



I. INTRODUCTION 

In studying a heat diffusion problem recently one of the authors (M.L.G.) encountered the following fractional 
_^' , integral 

^: i^rd.'^-^m±^ . (1) 



This integral does not appear in a succinct form in the more familiar tables of integrals such as Gradshteyn and 
Ryzhik [Q and Prudnikov et al although the latter gives a more general version of it in terms of 2^2 and li^i 
hypergeometric functions (see No. 2.8.7.6). So, a new method was devised for evaluating it, unbeknownst to the 
,—1 • authors at the time that the succinct form they eventually obtained was given on p. 139 of Apelblat's book [gj. 
J> ' Since this simple method turns out to be very effective for evaluating fractional integrals, we aim here to present the 
00 . technique with many new results to demonstrate its utility. In addition, we aim to extend the method in different 
CN ' directions to illustrate its utility and in so doing, we shall present many results not published previously. 
As we shall see, our technique can be applied to fractional integrals of the form 



^ — • r°^ F(x) 

''pi ,, with the aid of tables of LTs such as those by Oberhettinger and Badii Q, and by Prudnikov et al g]- Q and 
C^ • extended to many other classes of integrals, not of a fractional form. The technique can also be used in the numerical 
computation of integrals, since a slowly convergent integral can be transformed into a rapidly convergent one. 



II. THE METHOD 



Theorem: Given that the inverse LT for F{x) in Eq. (g) exists and is denoted by G(i), then the integral Ifj,{p) for p 
real and greater than zero can be written as 

/"OO 

/^(p) = r(l-Ai) / dtt^-^e-P'G{t) , (3) 



provided Re /x < 1 and G{t) is 0(i") as i — ^ 0+ where Re a > —Re /i while as t -^ 00 G{t) is 0{t^eP*) where 
Re /? < -Re fi. 

Proof: Since the inverse LT for F{x) exists, Eq. (||) can be rewritten as 

/'OD -1 /"OO 

IM=J dx—-—j dt expi-xt)G{t) . (4) 



Taking note of the conditions on G{t), we can interchange the order of the integrations whereupon we note that the 
x-integral is simply the LT of {x — p)~^Q{x — p), where Q{x) represents the Heaviside step-function. Now from p. 21 
of Ref. we have 

/•OO 

Clix-pyeix-p)} ^ dxe~^\x-p)-'' = t''-^e-P*r{l~fi) , (5) 

Jp 

where Re /i < 1 and Re a; > 0.[ The condition on /i was overlooked on p. 17 of Prudnikov et al ||^.] Introducing Eq. 
(g) into Eq. (Q) then yields Eq. (||), i.e. /^(p) = r(l — fi)Cp{t'^^^G{t)], where the subscript p denotes the variable 
with which the LT is taken. 

Before developing further results from Eq. (j^) we now consider some examples. Let us first consider the known 
fractional integral 

Up)- I dx l_J^ , (6) 



p 



which appears as No. 2.2.6.24 in 0. Thus, F{x) = x°'{x + a) '^ and its inverse LT is given in Ref. y| as 

/:-^{a;°(a: + a)--}^ ^ ° '^ iF,{iy;^a + u;-at) , (7) 

L [1/ — a) 

where Re(a — v) < 0, Ke x > {0, — Re a} for a, —v ^ 0, 1, 2, ..., Re a; > —Re a for a = 0, 1, 2, ... and Re a; > for 
-V == 0, 1, 2, .... Then using Eq. (|) we find that 

IM-T^ -J dtt^+-'-"-^e-P\F,{,,;,,-a;-at) . (8) 

T{iy - a) Jo 

The above integral is given on p. 510 of Ref. B and hence, Eq. (pi) becomes 

, , r(l — /i)r(/j, + u — a — 1) , , . , . 

^Mb) = ^ 2Fi{i^,fi + i^-a-l;iy-a;~a/p) = (9) 

It^ip) = (P + a)" V"^+" Bil~^l,^i + ly-a-l) 2^1(1 - ^,,iy;iy ^ a; (l + p/a)-^) , (10) 

where < Re(l — /x) < Re(i^ — fi). 

For the special situation where a = 0, we find that Eq.(9) reduces to 






^.(p)=^^73^ ;.^" ::..:-; , dD 



where we have utilised No. 7.3.1.1 from Ref. Q. For /i = 1/2 and p = 1, the above result reduces to No. 3.196.2 in 
Gradshteyn and Ryzhik ||l|. For the special case where ^ = we find that /^(p) — I3~^^{v, (3/{a +p)), where 7(1^, x) 
is the incomplete gamma function, while if /.t = 1 — v /2 and Re v <2, then 

V^r(^/2)r<'--^)/^e-^/^(-+rt r (UL\ T ( ^ M i^^^ 

The above analysis has demonstrated that our theorem can be used during the intermediate steps of an evaluation 
to obtain important results such as Eqs.(9) and (11). To obtain results not given the tables of integrals, first put 
V = 2n + v^ yu = — /i — n in (11), then multiply both sides of the original integral by (— l)"/3"/n!. After interchanging 
summations and integrations, we get 

Jo "^ — TpPpv — - p-M-irH 2^2 (^/i + l,/^-M- 1, 2'^-.-4^J > (-L3) 

where Re /i > — 1. This is equivalent to the interesting formula 

\^(X- xT-^-^e-^^^'^Mx = ^(^+^^^1;^-^^-^^ 2F2(M + 1, . - , - 1; ^, ^; -^). 



Similarly, we also find 



Jo "'' (t+p)'^-° •^°' \ t+p ~ 2=p-'-A'-ir(a+i)r(iy) ^ 



2F3(A* + l,j/-//-l;|,i^,a+l;-^j , (14) 

/„' x^(l - .x)''-A.-2g-/3.^^^ = B(j/ - M - 1, /i + 1) 2i^2(^, f + 1; f , ^; -/?) (15) 

3^3V3' 3 ' 3 '3'3'3'/' ^ ' 

j^;:pjy 2^31;^+ 3,t^+ 3,;^ + -3-,*^ + -3-,!^ + -3-,-2 j U'j 



r(3i/-M+4)r(iy+l) 



jl x-\l- xy-^ ^F^{^-^-v,!^^^-z^x^) ^ 
p7^T27r r("+^) ^■^+^-i(z)J,.-M-i(z). (18) 

There are many inverse LTs listed in Ref. [pj, which have the form x°'{x + a)"". For example, if we put a = —v in 
Eq. (|6|) and use our theorem, then for Re t/ > we find that, expressed as a Mellin transform, 

Jo [{x+a){x+b)y'^^ - 

il^,r-'B{s, 2^ - .) 2F1 (^ - f , ^ - ^; ^ + i; (f-f )2), 

(19) 

where we have used No. 2.15.3.2 from Ref. 0. One can come up with many interesting special cases, such as 

-dx = 4-^r(. + i/2)r(. i/2)_ ^^^^ 



/o ^(x2 + a2)(a;2+fe2)J r{2v){a + bf^-^ 

To consider an example that does not appear in the standard tables of integrals [fl , , H , let 






Then utilising our theorem we find that 



«^'^^;^ z"-'^"^^-'-"'-^(-n'-T-^) 



(22) 



which is valid for Re a + 2;^ > and for Im \a\ = according to p. 28 of Ref. [||. Utilising No. 3.38.1.17 of Ref. [gj, 
we find that the above integral yields 

3F2(^+^^±f^,^+^,i.;i.+ f,^+2±i;_^) , (23) 

where Re (/i + a + 2iy) > 1. 

For the case of a = 0, the above results yield the Mellin transform 

Jo f(x+l)2+a21-«2; - 



2 ' ^ ' 2^1 "^ ' 2 ' 



B(s, 2i^ - s) 2Fi{v -^,v+^^-v+\- -a2), (24) 



which is vahd for Re < Re s < 2Re i/, and 

poo x'^~'^^'^ J 

Jo {x+l)[{x+iy^+ai]-"-^ = 

(^r%TT^(VlT^-ir2Fi(.,i;. + l;i±4±iI). (25) 

Many other interesting results can be obtained in this way; several are given in appendix A. 

Now if we return to Eq. ( |24| ) and put u = k + u, then after using the duplication formula for the gamma function, 
we get 



r(/c+i/) roo , {x^+a'^)-^''+'''> _ 2^~^r(l-p) r{k+iy+fj./2-l/2) 

T(k + V + fl/2) Jp "'^ (x-p)M ~ („2_|_p2)fc + v + M/2-l/2 r(fc+i/+l/2) 



f dx 
Jv 

2F,[k + v+>^,l^:k + v+l-jJ^) . (26) 

Multiplying both sides of the above equation by (7)fca^'^/fc!, where (7)^ = r(fc + 7)/r(7) and then summing from 
fc = to cxD, we eventually arrive at 



M. g" ^ - r(2;y+A'~l) r(l-p) 

(a2+p2)i-+M/2-l/2 r(2zy) ^ 



2^i(^ + ^,7 + i^;^ + i;^) . (27) 

In obtaining Eq. ( P7|) we have used No. 6.7.1.7 from Prudnikov et al |8|. We may also multiply Eq. ( p4|) with u = k + v 
by t'^ /kir{k + v + 1/2) and sum to obtain the Mellin transform 

poo j.a-1 p(E. ii^- ,/ -U 1- ^-^f^'* — 

Jo [(2; + l)2+a2]'- 2-fU2' 2 ' '^ "•" 2 ' (a:+l)2 / " 

2-"r(s) 2i^i(jy - |, 1/ + i^; ,y + i; 62 „ ^2) (gg) 

where |a| < 1, |6| < 1, \b'^ - a^j < 1 and we have used No. 6.7.1.11 from Ref. ^. 
Before considering a numerical example, let us now consider the following integral 

1= f dy f 2J^i(a,fe;c;-c^yVp') , (29) 

Jo (1 - vY 

where Im w = 0. By making the change of variable, y — p/x, we find that Eq. IGh) becomes 

f°° x^~'^~^ 
^^P"J '^^ (x-pY 2Fi{a,b;c;-Lj/x^) , (30) 

whereupon we notice that the integral is now in the form of a fractional integral. Thus after applying our theorem 
we find for Re a > that 

/= El (°° dtt°'-^ e~P* 9f^ (a b-c °+i-^ a+2-M . -<^t^ \ ^ 

EM ,p,(a b - SL±1- c "+1"^ "+2-^ . i^\ ('o-|^ 

We can use the above result in conjunction with our theorem and specific inverse LTs to develop some interesting 
results. For example, if we put & = a+l,c=2a and a — ^ + 2a, then Eq. ( P9| ) can be written with the aid of our 
theorem as 

Jo "y {l-Vr 2^1 \"' " + -^' ^"' ?3-J - 22-4a^a-l/2r(2a+l) JO "^ ' ^ 

e-^* j2_i/2(V^i/2) = M£) 3F2 (a, a + f , a + ii±i; 2a, a + i; -^) , (32) 

where Re a > 1/2 and we have used No. 3.35.1.20 from Ref. §]. The above is basically the LT of ti'Jl_^^^{uj^/'^t/2). 
Hence, our theorem can be used to develop new LTs as well as Mellin transforms. 
To complete this section we now consider the following numerical example 



H^m-J^ ax ^^2+^2)1/2(2. _p)p ' 



I.ip) = dx ,,2 :^: ""^ „.„ , (33) 



where Im a ~ 0. Then by using our theorem we get 

/^(p) = -r(l - fi) p dt t^-ie-P* (^Ya{at) +{c + hi '^\ Mat)] , (34) 

where C denotes Euler's constant, <Rc /i < 1 and Yo{at) is the Neumann function. We can evaluate the above 
integral, but the final form is not very useful as we shall see. By using Nos. 2.12.8.4 and 2.13.6.3 from Ref. B the 
above integral yields after a little manipulation 

Up) - '-^^^^^^^^ 



2+a2)f/2 

r 



(ln(a/2) - C)F,_i { ^,./,.y^. ) + Q,^i ( (,.+:.) v. ) 



where Pu{x) and Qu{x) are respectively the Legendre functions of the first and second kind. As can be seen the above 
result is awkward because of the partial derivative. Since Jo{at) is bounded and Yo^at) can be written as the sum of 
a bounded term and 21n(at/2), Eq. ( p4| ) is simply the sum of an integral yielding the gamma function and another 
integral of the form 

/•CJO 

dt f"-^ e-P* \n t , (36) 

which is rapidly converging at the upper end. The integral can be handled at the lower end provided Ke fi > because 
then the algebraico-logarithmic singularity can be numerically integrated by following the procedure as described by 
Piessens et al on p. 47 of Ref. [|0|. 

III. NON-FRACTIONAL INTEGRALS 

In the previous section we remarked that our theorem could be extended to integrals not necessarily of a fractional 
form. To see this, note that we can write any integral as 

dx Fi{x)F2{x) . (37) 

We should add at this stage that any integral with other definite limits can be converted to one with the above limits. 
Assuming that the inverse LT of Fi{x) exists, we can write Eq. ( |37| ) as 

/ = J^ dx F2{x) /„°° dt e-* C^\F,) = /,°° dt Ct{F,)C^\F,) , (38) 

where we have assumed that it is valid to change the order of integration. Now one can see that if the inverse LT of 
-Fi(a;) has an exponential factor in it as, for example, the inverse LT of {x — p)~^Q{x — p) has, then Eq. ( pq ) simply 
becomes another LT. However, this is not necessary , although such a factor is particularly advantageous in view 
of (1) the availability of tables of LTs and (2) its amenability to numerical computation. As a consequence, in this 
section we shall concentrate on inverse LTs of Fi{x) that have some form of decaying exponential behaviour, e.g. the 
MacDonald function, K^{x). 

Let us begin by considering a simple example, which does not appear in the standard tables of integrals 0|, @, 01 
but is a special case of an example considered in the previous section: 

°° x-^+'jx-p)'^ 
p (a:2+a2)3/2 ' ^''^^ 

where the integral is defined for —1 < Re j/ < 1/2. Since C^^{x/{x^ + a^)^/^) — tJo{at), the above integral can be 
written as 

/ ^ /„°° dt t Jo{at)Ct{x-^{x - p)n ^ ^^-+11^' X 

f;^dtt-''+'/^e-P'/^Mat)K,+,/,{pt/2) . (40) 

Eq. (Eg) is still unknown but can be evaluated by expanding the Bessel function into a series. Hence, by utilising No. 
6.621.3 in Ref. §, we get 



r(t/ + l)r(l-2z^) /I 3 V Z-v a"^ 
/= , „ „,- ^ — 3F2 - -I/, 1- i^, -;1 - -,— r— ; 5- 



^l-2Iy 



r(2-^) 



(41) 



Eq. ( pl| ) reduces to the known result of (a^ +p^) ^^^ for i/ = 0. The result can also be checked by putting /i ~ —v, 
V = 3/2 and a = -jy - 1 into Eq. (|||). 

We can use the modification of our theorem to display similarities in apparently different integrals. For example, 
consider the following integrals: 



,2_„2N-l/2 



and 



/i = \^ dx x-^'e-P'"' [x + Va;2 - a^Y + {x - ^x^-a^f {x^ - a^) 



J^ dx x-^e-l'/'' {x^ - a^Y 



(42) 



(43) 



For /i we require He fj, > Ke v + I while for /2 , Re /i > 2 Re i^ + 1 and Ke ly > —I. To evaluate both integrals we also 
require Ci^x't'e-f'/'') = (i//3)'^"^'/V^_i(2v^). Then we have 



h = ^d'ii^^-^^l'^ 



dt i(^-i)/V^_i(2y^)i^,+i/2(ai) 



and 



h - ^'"^t-i)/2^^ '^ [ ^' t^/--'^-^J,-,{2^t)KM) 



(44) 



(45) 



Now we can see that both integrals are of the same type, which is given as No. 2.16.22.6 in Ref. [EJ, and find that 



h = 2^-ia'^-^-i 



aB{i^ + \,^-\) X 



2-Fri 



ti+u I 1 ti-v „ 1. 1 M M+l . P' 



4' 2 



4 ' 2 ' 2 ' 2 '4. 



^)-2/3B(H±il + |,H^ + i) 



2-^3 1 2 '^4' 2 ~'~4'2' 2 ' 2 'it? 



(46) 



which is valid for Re a > and Re /x > |Re(z^ + 1/2) |. Eq. (Esf) becomes 



2^~^a^"~^r(iy+l) 



2-F3 



{ t — I it 
\2 4' 2 



g r(;j/2-l/4) r(;j/2-t/-l/4) 



1. 1 £ i£+l. ^ 



4 ' 2 ' 2 ' 2 ' 4a2 



r(M) 

2/3 r(;j/2+l/4) r(p/2-iy+l/4) 



2-^3 



^3Ni + i ^ 



1. 3 M+l M 



4' 2 ' 2 '2 



■^ + -^' 4a2 



r(A'+i) 



/3-^ 



(47) 



the latter integral being valid for Re /i > | Re z^|+ Re z/ + 1/2 and Re a > 0. 

We should also point out that occasionally a fractional integral, which may not solved by using our theorem can 
be solved via the modified approach presented in this section. As an example, let us consider the following fractional 
integral 



/ = / dx {x- pYx'' [x^ -px + b) 
Ip 



-1/2 



(48) 



which is defined for — 1 < Re /i < 0. To evaluate the above integral we shall also require that | Im (6 — p^/4)^/'^|— Re 
p/2 < 0. In order to utilise the theorem given in the previous section we require the inverse LT of x'^(a::^ — px + 6)^^", 
which is not given in Refs. Q and Q. However, the inverse LT of {x'^—px + b)^^/'^ appears as No. 2.1.6.1 in Prudnikov 
et al |q] and hence, modification of our theorem yields 



/ = 



_ pM+1/2 r(p+i) 



2 



/;= dt i-(M+i/2) j„(t,/b-^wj^) K^+,/2ipt/2) 



CfT r(M + 1) r(-/i) 2^1 (i -m; 1; 1 - |I) 



(49) 



In obtaining Eq. ( p9| ) we have used No. 2.16.21.1 from Ref. [g[. For the special case of /i = —1/2, we find that 
/ = 2p^^K((l — Ab/p^Y'"^) where K(a;) is the complete elliptic integral of the first kind. Finally, more results based 
on the modification of our theorem appear in Appendix B. 



IV. FURTHER EXTENSIONS 

In this section we consider some extensions of our study of fractional integrals, the first dealing with integrals of 
the form of 

/•OO 

dx {x - py {x - piy^ F{x) , (50) 

p 

where p > pi, Re /ii > and < Re /^ < 1. Utilising the integral representation for the gamma function, we can 
write Eq. ( [5^ ) as 

-1 poo />oo />oc />oo 

' dx / die-^*G(t) / dz zi"-^ e-^'^-P^' dy y"''^ e-^''-P'^y , (51) 



r(/i)r(//i) Jp 7o 







where G{t) is the inverse LT of F{x). After a change of variable and evaluation of some of the integrals, one arrives 
at 

1=^^ Io°°dte-'^^Git) J^°°dyy^^-^e-(P-P^^y{y + tr-^ = r(l - m) x 

(p_p,)-(M+m)/2 J-rfig-(p+p,)t/2 G(i)t(A'+Ml)/2-l W^^_,,^^,^^^+^,_^y,[{p~p^)t) , (52) 

where we have used No. 3.383.8 from Gradshtcyn and Ryzhik |jl| and PVy,p(z) denotes the Whittaker function. This 
result can be analytically continued to Re /i < 1 for any value of /ii provided the integrals in Eqs. (|5y) and ( [s^ ) 
are defined. The condition on Re /x arises from the factor of r(l — fi) in Eq. (p2[). If we had introduced the integral 
representation for the gamma function in our theorem in Sec. 2, then we would have found that < Re ^ < 1 to 
obtain Eq. (0). However, we have shown that our theorem is valid for Re /z < 1 by using LT theory. This condition 
on /i arises due to the r(l — /i) factor in Eq. (m but there is no lower bound on fj,. Thus, while using the integral 
representation for the gamma function is initially restrictive, analytic continuation can be used to extend results such 
as Eqs. (0) and (p2|) beyond these restrictions. 

As an example of Eq. (|5^ ) in action, let us put F{x) — [x"^ + ax + 6)^^/^, /xi = 1 and /j, = 0. Then the original 
integral given by Eq. (pfl) can be evaluated by using No. 2.2.9.41 of Ref. [Q. The inverse LT of F{x) is given as No. 
2.1.6.1 in Ref. g. Hence from Eq. (|5^), we obtain 

I = j^dtt-^/^ e-(P+Pi+a)t/2j^(i^52 _ aV4)I^-i/2,o((p-Pi)i) = 

/ P-Pi l„ (p^+gp+b) 1/^ {pi +api +b) ^/^+fc+ap/2+pi (p+a/2) , r^-. 

y p'i+api+b (p-pi)(pi+a/2+(pJ+api+&)i/2) > \^'^) 

where (pl+api+b)^/'^ > 0, {p'^ +ap + b)^/'^ > 0, {p + a/2) > -{p'^ + ap-\-b)'^/'^ andRep> |Im(62 _aV4)^/^| ~Re a/2. 
For {-pI - api - b)^^^ > and {p + a/2) > {p^ + ap + b)^/'^, we find that 



arccos ( 1 + Pi+-Pi+b ] (54) 



(Pi+api+b) y (p-pi){^p^+ap+b+p+a/2) 

We can also develop another integral representation for Eq. I p(\ ) by using the Feynman integral [Q, ||l^ which 
allows us to replace the denominator of Eq. ( pOJ ) by 

^ „. ^ „ r(/i+Ul) f^ , t^l-1 (l-t)^-l , , 

ix-p)-^' (x-piy^ = ;\ i / dt- --^ ^ , ^ , 55) 

^ ' ^ ' r(Ai)r(Mi) io {{p-Pl)t + x-pY+^^^ ^ ' 

where Re /z, fXi > 0. In Appendix C we establish the above result from its more general version, which appears as No. 
2.2.6.1 in Ref. 0| and then use these integrals to develop several interesting results. Thus, Eq. ( |5(3| ) can be written as 

I = ^f^Wl^ r dt --''m /o ds s'^^-Hl sr-^ /- dy (,^(,4:;)..., . (56) 

The final integral in Eq. (|56| ) can also be evaluated by using No. 3.383.4 in Gradshteyn and Ryzhik Q. So, Eq. (|5(| ) 
becomes 

/ = ^f^^^ /o" dt e-P*G(i)i(^+^i^2)/2 Jl ^g sPl-l(l _ s)/^-i((p - pi)s)-(''+A'i)/2 X 

e^^-^^^^*/'W^-(,.+M0/2,(i-M-Mi)/2((p-Pi)5i) • (57) 



If in Eqs. ( |5^ ) and ( p7| ) we put G{t) = Jo{zt), i.e. F{x) = [x^ + z^) ^/^, then by multiplying both equations by 
zJo{zy) and integrating over z from zero to infinity, we obtain the interesting result of 

r(/.i) e(P-P^)/2 iy(^_^,)/2,(Mi+M-i)/2 (P - Pi) , (58) 

where < Re /x < 1 and Re /ii > 0. 

For the special case oi ji = jii the integral given by Eq. (BG) reduces to 

/;° dx (x^ -{p + p,)x + ppi)-^f (x) = ^ (^i^-^g-v. /„°°die-(^+P^)*/2G(t) X 

i^-i/2if^_^/2((p-Pi)i/2) , (59) 

where we have used Eq. (52) and the fact that lVo,^(z) — ^J z / t: K ^^^{z / 2) according to No. 9.235.2 of Ref. ||l|. Eq. ( |59| ) 
can be checked by putting F{x) = x~^ whereupon the integral on the Ihs after a change of variable becomes 

i={p-p,)-^p^-^-^^^dss-^(i+^y\i+s)-'^ = 

{p^p^)-^^p^-^^-B{l-^l,u + 2^l-l)2F^(^l^-^l■^l + y■,-^) . (60) 

In obtaining the above integral we have used No. 3.197.5 from Ref. [nl. It is valid for Re /i < 1, |arg(p/(p — Pi))| < ti" 
and Re i^ > 1 — 2Re /i. One can show that the rhs of Eq. (|59| ) gives the same result by noting that Lf^{x) is t'^~'^ /T{v) 
and utilising Nos. 6.621.3 and 9.131.1 from Ref. ||l|]. Furthermore, if we put pi — ~p in Eq. (|5^), then we obtain 

j;^dx {x^-p^)-^F{x) ^m^J^^dt {t/2pr-^/' G{t)K,-i/2{pt) , (61) 

which can also be found by following the procedure outlined in the previous section where one would eventually 
evaluate the LT of (x^ — p^)'^ ■ Using the latter approach, however, means that the restriction on /i as given under 
Eq. (^) can be analytically continued to Re ^ < 1. 

A pleasant example that can be derived from Eq. (p9|) is 



-1 

(2a-l)^ 



B{l+u,-2u) f a-c p/1 -,.-, .,. ^ a(a-l)+2e(l-c) N 
^T^ 12^^ 2^^1(2: i:-L-J^, 4 



6- 



' P('l 1-1 ,,. 4 b(b-l)+2c(l-c) fnr)\ 



26-1 2^ l\2^ ^ ' (2 

where, to ensure convergence, — 1 < Re v <\. In obtaining this we have utilised various results given in Refs. |0] and 

To set the scene for the material to follow let us examine Eq. (|61|) more closely by analysing the behaviour at the 
lower and upper ends of the integrals on both sides. Clearly for the integral on the Ihs to be defined, F{x) must be 
0{{x—pf) as X -^- p^ where Re j3 > Re /i— 1. This condition has been implicitly assumed throughout in discussing all 
fractional integrals. At the upper end, F{x) must be 0(x"), where Re a < 2Re fi— 1. If we were to replace /x by /i — 71, 
where n is a positive integer or equal to zero, then the Ihs of Eq. (^ would only be defined if Re a < 2Re ^ — 2n — 1. 
Thus, if we wanted to sum over all values of n from zero to infinity, then to guarantee convergence F{x) would have 
to possess some exponential decaying factor such as exp(— ex) as a; — > 00, where Re e > 0. However, this may be too 
restrictive for many integrals may be defined after the summation process has been effected. 

The singular behaviour is reciprocated when examining the rhs of Eq. (ply. Thus, if we were to replace n hy fj, — n 
in the rhs of this equation, then for the integral to be defined for jjl — n— 1/2 < we see that G(i) must be 0{t°') 
as i — + 0+, where Re a > 2n — 2Re /i. Again, we can overcome this singular behaviour by requiring G{t) to have an 
exponential factor. In this case it would be exp(— 7/i) but it would also be too restrictive on G{t). 

By taking note of the above, we can introduce exponential factors to both sides of Eq. ( |6l| ) so that it becomes 

lim,_o+ r(l - ^)-i J^ dx (x^ - p2)-pg-e(x-p)^(^) ^ lim,^o+ /o°° dt e~''\t/2pY-^/^ x 

7r-i/2 G(i)i^,,_i/2(pi) . (63) 

This enables us to replace ^ by /i — n and then to sum from n = to 00. As a consequence of absolute convergence, 
we can interchange the order of the summations and integrations. Finally, we can take the limit as e -^ 0+, which 



will yield integrals with less restrictions on the behaviour of F{x) and G{t) than those mentioned previously. Thus 
by putting /i = ^ — n and multiplying both sides by (— l)"/3^"~^/2^"~'^n! of Eq. (|6^), we obtain 



l™e^o+ EZo 2^^^^^n) ir d- {-' -p'r^^-e-^^-'^Fix) ^\im,^o^j:Zo¥^ x 



Since the integrals are now absolutely convergent, the order of the sums and integrals can be interchanged to yield 



/p°° dx (x^ - p'r/^J^iPy^^^^^^)F{x) = /o°° dy -^ J^{Pv)nV¥+P^) 



1 pM-i/2^M ;- dt {t^ + f3^)-^/^-y^G{t)K^+,;2{pV^TW) , (65) 

where we have taken the limit after evaluating the sums and have also replaced /i by —/i, so that the above result 
is valid for Re /i > — 1. In addition, we have utilised the multiplication theorem for the MacDonald function, which 
appears as No. 14 on p. 112 of Ref. B. 

Let us test the result given above by putting F{x) — x^^'^^^ . Then the Ihs of Eq. ( pSJ ) becomes No. 6.565.4 in 
Gradshteyn and Ryzhik ||I|, which is 

This result is valid for — 1 < Re /i < Re (2;/ + 3/2) with both /3 and p greater than zero. According to No. 2.1.1.1 in 
Ref. [pi, the inverse LT of x"^""^ is t^''/r(2j/ + 1), which when introduced into the rhs of Eq. ( |6^ ) yields the integral 
given by No. 6.596.3 of Gradshteyn and Ryzhik. This, in turn, yields the rhs of Eq. ( pq ) after application of the 
duplication formula for the gamma function and is valid for Re v > —1/2. Thus the integral emanating from the rhs 
of Eq. (|6^) represents the analytic continuation of the integral on the Ihs of Eq. ( pS] ) . 

We now consider an example, which does not appear in the standard tables of integrals. This is 

/ = /o°° dt {e + /32)-M/2-i/4i4^ ^p^ (^^ + i; 2a. + i, 2m + 1; ^) K^+,^,{pV^TW) ■ (67) 

Utilising No. 2.1.7.20 of Ref. |q| and Eq. (|65|) above, we can transform Eq. (pT]) into 



/ = 



f j^^ r dy y-^y' + p')-^-'^' MM = ^f r(M + 1) x 



2'- r(4M+i) J §p\j^ §p\ (68) 



where Re /i > -1/2, /3 > and we have used No. 2.12.4.30 from Ref. ||]. 

An interesting application can be found by putting y^ = — 1/2 in Eq. (pa), which yields 



^^dy^2^^F[,/fTW)=p-'S^dtG{t)K,(p,/WTW) . (69) 



If we put G{t) — Ki{p{t'^ + p^Y/'^)/{t^ +p^)^/^, then using the result in the appendix of Ref. |13| we find that the 
integral on the rhs of Eq. (||) yields 'KKo{2pl3)/2(3p^. We can identify F((y2 +p2)i/2) ^y noting from No. 6.726.4 of 
Gradshteyn and Ryzhik p| that 



37r /-oo 

2/3 



■ r dy ^.^:^)^,. ^1/2 {N^iv^+P')) = 2k^ K,{2pP) . (70) 

Hence, we obtain the following LT which does not appear in Refs. Q and g]: 

/o°° dt e-y' {e + (i^f/^K^{p^WTW) = y5p-2^l/2(^/3/3y) . (71) 

We should mention that the approach used to derive Eq. (p3) from Eq. (pi]) can also be applied to Eq. (p9), 
whereupon one obtains 

/;" dx {x - pY/\x - p,Y/^J,{fi^{x-p){x-p,))F{x) = ^'^•'-1%""'' X 

/~ dt e-(P+Pi)*/2G(i)(t2 + /3')-'^/'-^/"i^^+i/2 (^^v/t^+^) ■ (72) 



In the above result we have put /i — —fi, so that it is vahd for Re ^ > —1. Eq. ( Jl2( j can be verified by multiplying 
both sides by /3J^(/3y) and then integrating over /3 from zero to oo. This gives the correct result of F{x) as the LT 
of G{t) after utilising No. 6.596.7 from Gradshteyn and Ryzhik 
making an appropriate change of variable, one can express Eq. 



. In addition, putting pi = —p yields Eq. (|65|). By 
72h in the more convenient form of 



(,2+^2)-,/2-l/4^^^^^^(^fn^^ . (73) 



In Eq. (jT^) one can replace a/2 by a and a/2 + pi by p, which we do henceforth. Multiplying both sides of Eq. ^ 
by P"^ and integrating over /3 from zero to cxd, we get 



r{{f,-^+l)/2) a('--'-)/^ ^ 



C ds f^ F (V«2 + s^+pj- ^-^ 

2(M-.)/2 J- ^^ i(.-p)/2 e-P*G(t)i^(^-.)/2 («t) , (74) 

where Re {p + u) > —1 and Re i^ < 1/2. 

We now apply the above procedure to obtain an apparently new result. Consider the following fractional integral 



1 = J^dt e-P\t - a)^-i/2 (^2 ^ p-^yt^^K^ipy/WTW) ■ (75) 

Comparing the above with Eq. (^ we see that G(t) = (i — a)^^^/^ and hence, from No. 2.1.2.10 of Ref. [||, 
F{x) = r(^ + l/2)a;-(^+i/2) exp(-aa;). Thus, Eq. (||) yields 



/ = 



-" r(M+l/2) TOO , +1/2 . _ ^^,^ X e-°v'^^+^ 



- r(/. + l/2) ^l/2^-p(v/a2+/32 + „) 

2 P^ (Va2+/3Ha)f'-i/2^a2+/32 

where Re /i > —1/2 and we have used No. 2.2.10.13 from Ref. 0. 
Now consider the following integral where a > and 6 > 0, 



+s^ 



(76) 



i = ^rdy^-^^^{y'^b-)'^-'^--'/' . (77) 

If we make the change of variable, x = [y^ — 6^)^/^, then the integral is in a form given by the Ihs of Eq. (ItI) where 
F{x) = x~'^l'^ exp(— a2/4x). By utilising the fact that the inverse LT of F(x) is given as No. 2.2.2.3 in Prudnikov et 
al 1^, we find that Eq. ( ^ ) becomes 

2r((^-i.)/2+2) 1-^1 (^'^"Ai+ 2'^-+^'^8hJ ' y'^) 

where |Re (p - v)/2\ < Re (i^ - ^ + 3)/2 and we have used No. 2.16.18.2 from Ref. 0. 

Let us return to the theorem given in Sec. 2 and examine the denominator of the fractional integral as given by Eq. 
(0). We see that an inversion in the power /i, which we aim to exploit in the remainder of this paper, occurs when we 
evaluate the integral by using Eq. (ph. So far, we have been concerned with situations where the integrals have been 
defined but what happens if we replace ^ by /i — n in our theorem and let n range from zero to oo. Clearly, it can be 
seen that the rhs of Eq. (|^) can be become divergent depending, of course, on the behaviour of G{t). To overcome 
this potential problem, we can use the work of Lighthill |1J] and Ninham [^, who give interpretations to the class of 
divergent integral encountered when we consider a range of values for n after substituting /z by /u — n in our theorem. 

If we put /i = /i — n in our theorem, then after multiplying both sides by (— l)"/?"/n! where Re /3 > 0, we obtain 

;nte^ r d^ Ij5S- - '-^ r dt i— e-^'G(0 . (79) 

In terms of generalised functions the integral on the rhs of Eq. (|79|) can be written as 

/ = /^ dt t^'-"-^e-P*H{t)G{t) , (80) 
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where H{t) is the Heaviside step-function. Lighthill jlj] shows that if t'^~^G{t) is a good function, then the above 
integral yields the ordinary formula for repeated integration by parts if all the infinite contributions from the lower 
limit of the rhs of Eq. ( |79| ) are omitted. As a consequence within the context of generalised functions, only the finite 
part to a divergent integral need be kept. This finite part is referred to as the Hadamard part by Ninham [[l5[ , 
since Hadamard showed that the finite part obeys many of the ordinary rules of integration. Ninham's approach to 
divergent integrals of the form given by Eq. (gO|) is different from Lighthill's in that he is able to express such integrals 
in terms of a sum of function evaluations over an arbitrary partition of the integration interval, thereby making a 
connection with the concept of a Riemann integral as a limit sum. Ultimately, his approach coincides with Lighthill's 
but for it to be applicable to integrals such as Eq. (|8^), one needs to make a change of variable such as y = p/{t + a) in 
this integral. Then the singularity is shifted to y = l/a. As a numerical example, Ninham chooses the beta function 
integral 

/(a,/3) ^ J^ dt r(l - tf ^ "'°+ WD . (81) 

According to his prescription, this integral should yield a result for all a and /3 and for convenience he puts a = —3/2 
and P = —3/2. The Hadamard contribution is the final expression on the rhs of Eq. ( pl| ) and yields zero for these 
values of a and (3. Ninham shows that by treating the integral as a sum over each and every partition of the interval 
[0,1], denoted by i?/, and an error term associated with the partition, denoted by Ef , the integral virtually yields 
zero. A slight error arises from the rounding-off of Rf and from the truncation of the asymptotic series for Ef . 

From the above, we can see that provided it is defined, the Ihs of Eq. ( |7g| ) represents the Hadamard contribution 
for the integral on the rhs. If we sum both sides of Eq. ( p9| ) over n from zero to oo, then after interchanging the order 
of the summations and integrations we get 

/o°"dii-^-ie-P*-''/*G(t)-2r^/2/J^dyy^+ij^(2V^2;)F(y2+p) , (82) 

where we have put x — y"^ + p and Re /i > —1. For the special case oi fi — —1/2, Eq. (p3) reduces to 

/^ dt i-i/2e-P*-/5V4t G{t) = ^ /„°° dy F{y' + p) cos{(3y) , (83) 

where we have put /? — /3^/4. Note that in the above results our summation procedure has transformed the singularity 
at t = from either a pole or a branch point, depending upon the value of /i, into an essential singularity, which has 
a less harmful effect on the evaluation of integrals. We shall see shortly that the summation of divergent integrals can 
also lead to a shifting of the singular point away from the integration interval. This technique of summing divergent 
integrals to obtain a convergent result is the essence of the renormalisation or regularisation technique used so often 
in theoretical physics. In addition, if we had multiplied by f3 instead of — /?, then (3 would have to be replaced by — /? 
on the Ihs of Eq. ( p2|) whilst the Bessel function on the rhs would have to be replaced by a modified Bessel function. 
To convince the reader of the validity of the above results, we now consider some examples. Let us put G{t) 
equal to unity, so that the Ihs of Eq. (g) yields 2{l3/p)~t'^^Kf,{2y^) by No. 3.471.9 of Ref. §. Since G{t) = 1, 
F(y'^ +p) = [y"^ +p)^^ and hence by using No. 6.565.4 of the same reference, the rhs of Eq. ( [8^ ) can be shown to give 
the same result as the Ihs. On the other hand, if we let F{x) = x~'^ where Re ;/ > 0, then the Ihs of Eq. ( |8^ ) yields 

/„°° dt i-3/2,-pt-,v47r(.) = ^^;^:^:;;;;^^ K^-i^i^^) . m 

To obtain this result we have used No. 2.1.1.1 from Ref. Q. Introducing (y^ -\- p)^'^ into the rhs of Eq. |8^ gives the 
answer above after utilising No. 8.432.5 from Gradshteyn and Ryzhik [|l[. 
A results obtained from Eqs. ( p2|) and (|8^), which appears to be new, is 



[Re p. Re /S^ > 0, Re i' > -1/2] , (85) 



To study a situation where the summation of divergent integrals produces a singularity shifted away from the 
integration interval we multiply Eq. ( |79| ) byr(r7, — /i+1), replace /x by —/x and then sum over n from zero to oo. Then 
after interchanging the summations and integrations we get 

/p°° dx [x - p)^F(a;)e-'5(--p) = r(l + ^) ^^ dt e-P'G{t)t-t'-^ i^o(l + M; -P/t) , (86) 
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where Re ^ > —1. Introducing No. 7.3.1.1 from Prudnikov et al M into Eq. (|8q) yields 



Git) 



(87) 



/;° dx [x - pYF{x)e~P^^'P) = r(i + ^,) S^ dt ^ 

Hence, we can see that the singularity has been shifted to — /3 and is now outside the region of integration. Eq. 
can be verified by letting F{x) = x~^ or G{t) = 1. To show that the Ihs is indeed equal to the rhs, put y = x — p to 
get 

which is just the rhs of Eq. (|8^). By multiplying Eq. ( |79| ) by r(n + v + I), we can generalise Eq. (86) to 

J;;°dyy^^F{y+p),F,{,,+ l■^^+l;-py) = ^{^^ + l)J;;°dtG{t)'^^^^ , (89) 

where Re /i > —1. In addition, if we put fi = n + fi and multiply both sides of Eq. (|87| ) by (— l)"7^"+^/2^"+^n!r(n + 
/i + 1), then after interchanging the summations and integrations we obtain 



pt--,-'/4,(t + l3) 



j;rdss^+\j,ijs)e-^^'Fis^+p) = ^S^dtG{t) ^^^9^ 



(90) 



We can still go further with Eq. d79). By dividing both sides by r(7i -\- v -\-l) and then carrying out our summation 
procedure we find 

j^ dx [x - pYF{x) 0F2 i^i+l,iy + l; -P{x - p)) - /?-''/2 r(/i + 1) T{v + 1) x 

S^ dt iM--/2-i e-p/t G(t-i) J, (2V^) , (91) 

where Re /i > —1. By putting F{x) — x^^ and using No. 2.12.9.14 from Ref. [|| with the above result we get 

/o" dy (^ oi^2 {l + ^^^ + v■ -py) = T{v + l)V(p + 1) [r'' r^T^ ^ 

(92) 



^F2{1- ^l^ + u- ^l■|3p)+p^^ ^^ oF■2{^i+l,v+l■pp)\ , 



r(i/+i 

where Re (2^ - j/) < 3/2, Re p > and V^ > 0. 

To complete this section we utilise the preceding material to obtain the following general result for Re /i > — 1 

!!^ dx {x-pYF{x) pFq+i {ai, ...,ap; fi + l,ji, ...,jg; -/3{x - p)) =r{fi + l) x 

/o"dit-^-ie-P*G(i)j,Fg(ai,...,ap;7i,...,7,;-f) . (93) 

This result represents the generalisation of the theorem presented in Sec. 2 and is not always restricted to Re /3 > 0. 
Using the tabulations of hypergeometric functions given in Ref. [pj , we obtain the following results from the above 
equation 

J^ dy y3/2-& F(y +p)Jfc_i {(3^) J^-b {P^) = p-' ^ J^ dt t^-^ e-^* x 

G{t)(l- ii^i(|;6;-f 



and 



[Re 6 > 0] , 

/;° dx F{x)Jl^,,, [VM^^) = r dt t-' e-^^-PI^^ Git) /,_i/2 (-1) 

[Re a > -1/2] . 

Putting b = 3/2 in Eq. (M) eventually yields the result given by Eq. (B3h. 



(94) 



(95) 
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V. CONCLUSION 

In this paper we have presented a theorem/technique derived from Laplace transform theory that facihtates the 
evaluation of Weyl fractional integrals by transforming them into known integrals. As a consequence, we were able 
to present many new fractional integrals not previously evaluated in the standard tables of integrals, Refs. |jl|- P|, m 
and |q]. Some of these results appear in Sec. 2 where we introduced and proved the theorem. Further examples are 
presented in Appendix A. In addition, the technique is particularly useful for the numerical evaluation of fractional 
integrals since it is able to transform slowly converging integrals into rapidly converging ones. 

We showed in Sec. 3 how the theorem could be extended to integrals not of a fractional form but concentrated upon 
integrals which yielded decaying exponential-like behaviour after transformation. As a result, we were able to evaluate 
many new integrals, which are listed in Appendix B. In Sec. 4 we considered several extensions of our theorem. During 
the course of our study we found that by utilising Feynman's integral from quantum electrodynamics together with 
the techniques in this paper we could derive an interesting class of results, anumber of which are given in Appendix 
C. Finally, by using the divergent integral theory of Lighthill WM and Ninham |la, we were able to show that our 
theorem was only a special case of the result given by Eq. (BSh . 
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VI. APPENDIX A 

We begin this appendix by utilising the theorem in Sec. 2 to evaluate a few more fractional integrals. In listing 
these results we also indicate which references were used, so that should any errors appear, the reader may trace their 
origin. 
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J^ dx iV^ + arix - pY 






2F1 (-/. _ I - 1, -^ - £±1; -/, - J.; 1 _ i^) 



[Re fi > -1, Re i^ < 0, Re ^ < -Re i^/2 - 1, |arg a\ < 37r/4] 
[No. 2.1.7.23 Ref. [|, No. 2.11.3.2 Ref. |]], 



(96) 



t^-M) 



f °° ^„ ( {V^^+^±xr \ , _ w _ g" r(M+i) r(Tt 

[Re /I > — 1; =FRe i^ > Re /x, p > |Ini a 
[No. 2.1.9.15 Ref. [|, No. 2.12.8.4 Ref. |]], 



r(^ + l)r(. - ^) ^(j^ 2F1 (.-/.,.+ i; i - A^; H^) 
[|Rei^| < 1, -Re^ > |Re v\, Re (a+p) > 0] 

[No. 2.1.9.24 Ref. [|, No. 2.16.6.3 Ref. |]], 



(97) 



/p°° c^a; (x - pY 



,1/3 _ r(A.+i) 



ri^7rTT^r(^+f)r(M- 



nt2 \-fl- ^,-fi- o,-^TT 



1 . paM , S-'p' 



3' 27 j "^ 2aV2 



r(-i)r(-A.- 



8 ±._P2l\ a_ 3V1!^ r 



0^2 (m 

[ 

[No. 2.2.1.4 Ref. §, No. 2.16.8.13 Ref. §] 



3' 3' 27 I ' 2a''/^ 

arg a\ < tt/S, Re /i > — 1] 



(i)r(-M-|)o^(M+i,|;-#) 



(99) 



r rf^ ^^ e-"^"''' = %^ k''+i/2r(-(^ + 1/2)) X 



P A. 3 ,, I 3. a^p 
1-^2 1^2' 2'/^+ 2' ^ 



2v^ 
2^^'r(^+l/2) 



ii^2 (-A^; 1 - M, i - /^; T^ 



[Re a^ > 0, -1 < Re Ai < 0] 
[No. 2.2.1.16 Ref. [|, No. 2.8.5.15 Ref. |]], 



(100) 



/7 dx {x - pY e''^-^erfc(ax) = H^tJ^;-!}^ e^^^^^^ ^^ (V2ap) 
[-1 < Re //< 0, Re a^ > 0] 

[No. 3.7.3.2 Ref. §, No. 2.3.15.3 Ref. §], 

r dx (2_^ (C.2 + ^2)2a-c+l/2 ^^^ („, ^ + 1/2; c; -u^ / x^) = ^,lttlU >< 



r(2c-2a-p-2) 



P 



2c-2q-^i-2 — 2 



(2c-2a-l) 

i^i(c - a - ^/2 - 1, c - a - ^/2 - 1/2; c; -uS^jp^) 



[Re c > Max[Re(Ai/2 + a) + 1, Re a + 1/2], p > |Ini w[ 
[No. 3.35.1.22 Ref. [|, No. 2.12.8.4 Ref. §], 



(101) 



(102) 
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f dx 
Jp 



{x-pY 



2n _ r(M+i) 



,^i(a,a + l/2;c;-^/x^) = ii^ x 



r(2a-^i-l) 
pia-M-1 2 



Tji / 2a — u— 1 2a^ 



2 ' ' p2 

[Re a > 0, p > |Im ^|, Re (2a - ^ - 1) > 0, Re /^ > -1] 
[No. 3.35.1.17 Ref. [|, No. 2.12.8.4 Ref. §[], 



(103) 



.,/^ r(7+l/2) ^ p (^^ ^ 7+1/2 7+»-l/2 .g ^.-c^n"") 



[r + n<s+l, Re 7 > 1/2, a, ^ 0, -1, -2, ...] (104) 

[No. 3.36.1.11 Ref. [|, No. 21.26 Ref. g], 

We also mentioned in Sec. 2 that our theorem can be utihsed in the construction of fractional integrals when it 
cannot be used to solve them. As an example we cited Eq. (O), whose more general form, given by Eq.(9), was 
evaluated by using our theorem. We were able to use this result to develop some fractional integrals. We should add, 
however, that more fractional integrals can be developed from this result, which we present here. Thus, by putting i' 
equal to n + z/ in Eq. ([ll|), multiplying both sides by (3"/nl and summing from n = to oo, we get after interchanging 
the order of the summation and integration 



where the same conditions apply to /i and v as in Eq.(9). Similarly, one finds 



rdx^^ 

Jp ix 



jx+a) 



{x-pY 



exp 



{'^)- %7il%tx:n ^FAi-,-,.-,P) , 



(105) 



(106) 



f dx 

Jn 



{x+ay 



{x—p)f^ 



Jl 



(m+i^)/2-1 ^+^ 



V^\ _ 2''-lp{,^ + :y}/4-'l/2 r(i-^) 



r((p+iy-i)/2) p / M+i^-i . 1^ v+i 

r((i/+l)/2) 1-^2 I 2 ' 2' 2 ' (p+a 



(p+ajf + ^-i 

/3 



TW 



(107) 



2^~^--T{{l-p)/2) 



r(i./2) 



r((l-M)/2) ^CL /l _ i£ 1 _ ii. ii idi- -R\ 

2-* 2v9 9' 979' 9 7 r^/ ' 



2 ' 2' 



Jp "-^ (2;-p)i/2(^+a)3/2 - V2 (p+a) '^'^/2 VSlP+a)^ '^('^-l)/2 V2(p+a) 



(108) 
(109) 



and 



Jp "-^ (k-p)a' 0-^3 1 2 ' 2 i*^' 4(a;+a)2 






r(i^+M~ 



Vtt (p+ 



arH-(5jT.)+^-^- 



/3 



2\/p+a 



Jp "^ (a:-p)f 






r(i-M) 



/a;+a 



'■I'+p-l 



/x-\-a 



T{v)V(u+ti-l) 0-^3 



Tii ii+i 
\^2' 2 ■ 



i^ + M - 1; 



/s-' 



64(p+a) 



(110) 



We also stated in Sec. 2 that more results could be determined from the integral given in Eq. (E] 
tabulated results for ^F2 hypergeometric functions in Ref. [B. For example, if we put a = —2 in Eq. ( pl[ 

r ^^ ^^^-P)^ = P'''-'" B{1 - A^, M + 2^ - 3) X 



(111) 

by using the 
, then we get 



2-f 1 



(. + V, 



M 1 . ,, 1 . a 
2 " -^'^^ 2'^P 



(t/+(M-3)/2)(t/+p/2-l)a^ 
(^-1)('^-1/2)P^ 



2F1 Z/ 



M-1 



(112) 
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If we put a ~ 1 and fi = 1/2 — v, then Eq. ( pT| ) yields 

r d- ^9^$^ = P-"^-^/^ S(- + 1/2, - + 1/2) (i±Vi+ZZZ) " 



17 / 1 ,1 l + \/l+a^V 



(113) 



Eq. (U13) is based on using No. 7.4.1.13 in Prudnikov et al M, which has another representation given immediately 
above it in the same table. Hence, the above result can be expressed in a different form. On the other hand, if we 
put V = 1 and a = 2, then we find 



°° ^^ x-\x^+a^)-' _ 



L ^^ 



{x—p)^ 



' (M^+3M+2)a^ 



iFl 



I 2 ' 2 



1; 



- 1 



(114) 



To complete this appendix we now give further results that can be derived as a result of the evaluation of Eq. (|2S|), 
whose answer is expressed as a 4i^3 hypergeometric function in Eq. (pi]). Thus, using the results in Sec. 7.5.1 of Rcf. 
fil, one obtains 



( 



y^\ _ -i Lj-'^/^p{a-fi)T{a) 



Fi [2a— 1, a — 1; a — ^; 



2(2a-l)(a-l)r(Q+l-;i) 



'"'" - .2F1 (2a - 1, a - 1; a - ^i■ -^ 



(115) 



and 



/o^ rfy igr 2^1 (f , 1 - M + f ; ^ ~ m; -1^) - "l?ilT!^;r 

2F1 (|, ii±l - ^; a + 1 - ^; z) 2^1 (|, S±l ; a + 1 _ ^; z) , 



/„^d2/laarctan(^) 



- 3t^'^^ ^-3/4 
~ 4p -^ 



(l + zi/4)2in(i + ^i/4) 



(l-zi/4)2ln(l-zi/4) 



f^ ln(l + VI) + ^ (V^ - 1) arctan(zV4) 



2p^/z 



2pz3/ 



In Eq. ( 116 ) z is a solution of z^ = —4uj{z — l)/p^ while in Eq. (117) z — exp(i7r)tj/p^. 



(116) 



(117) 



VII. APPENDIX B 



In this appendix we develop further results based on the material presented in Sec. 3. We shall concentrate on 
integrals where part of the integrand has an inverse LT with some form of exponential behaviour. As in the previous 
appendix the list of results presented here is by no means exhaustive. 

We begin with integrals containing the factor of x'^ exp{—a/x) in their integrands. Since the LT of this factor is 
given as No. 2.2.2.1 in Ref. D, we find 



/o°° dx x''e-^'^F{x) = 2a(''+i)/2 J^ dt t-(-+i)/2 K,+i{2^t) G{t) , 



(118) 



where G{t) is the inverse LT of F{x) and Re a > 0. We now present some new results arising from Eq. (118). As 
before we list the references used in obtaining these results, so that the reader may verify their correctness. Thus we 
obtain 

j^ dx :,-e-/--P-"" = [-P 0^^+3/2 F {-v - |) 0^2 [%v + |; -ap") + 
4a''+i V{~v - 1) 0F2 (i, i/ + 2; -ay") + -f^ V{v + 1) x 



r(i^ + 



Fo (-V 



0-f^2 



-ap2)] 



[Re p2 > 0, Re V^ > 0] 

[No. 2.2.1.9 Ref. [|, No. 2.16.8.13 Ref. @], 

By using No. 2.16.5.4 from Ref.[2], we get the new Laplace transform 



(119) 
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^^ -i y 2 ' 2 ' 2 



Ai-A+1 , -, ^-A+2 . g^ 
2' 2 '" ' 2 ' 4 



££ r(A-..)r(^./2)a^-^ 

4 [ r((M-^+l)/2) ^ 

r((A+i)/2) 



2^3 



Fs 1- 



A+1 ££+1.3 ti-X+2 3+£i-A . g^ 



2 ' 2' 



r(p-A)r((p-t.-l)/2) 



P (_££ A. j_ 1- f. , -. - I -. f. . ^ , — yr- -- . -,~ ^v- . -„ -, y 

2-^3 \^ 2' 2' 2' 2 ' 2 ' 4 J 2f'-^-ir((/^+l)/2) 



-/j+A 2+A--£t . 
2 ' 2 ' " 



4 y ' 2f'-A r(/i/2+l) 

gr(p-A+l)r((A+l)/2) 



3 2+A-^ 3-£t+A . g^ 
'2' 2 ' 2 ' 4 



2F, f i^ ^±i- 
2^312' 2 ' 

[Re /i > 0, Im a = 0, Re Va > 0, Re (^ + A) > |Re (/i + A|)] 

[No. 2.1.9.15 Ref. [|, No. 2.16.22.2 Ref. §]. 

Integrals containing a;'^(a; — a)"^ where a > can be transformed as follows 

C dx x-{x - ay Fix) = ^^-+^y^^"' /- dt t—^/\~-'/^K,+„MmG{t) 

where Re v > —\. Some examples arising from this result are 

r dx x^{x - ay '-^^>"^ = 2--l-a-^^-l% ^'^^^gL^^,'^-^ x 

/"s 5 1 

\^4' 4' 4 



4^3 



3_.3 3_££ 5_i£._b^ 
"^i 4 '^J 2' 4 2 ' 4 2' g2 



[&>0, -1< Rejy < 1/4] 



[No. 2.1.8.25 Ref. [|, No. 2.16.18.1 Ref. 



(120) 



(121) 



(122) 



\^dx x^(x- a)''exp(-6a;i/2) = £(ig 



r(t^+i) 



/ — 2 

yTT g 



22,^ + 2 



i:^r(. + i)r(2.- 



ii^2(i^ + l;2i. + 2,i;2f 
i^2(-.;-2.,-2.-i;^)+™-+3/2,ngl-^ X 

iF2(z.+ |;2^+|,|;^)^ 
[Re62 > 0, Rei^ > -1] 

[No. 2.2.1.9 Ref. [|, No. 2.16.9.3 Ref. |]], 
Integrals of the form 

^-r^^v^ {[x^^^^^^y ^{x-.fx^^^y\ 

can be written as 

I ^ 2a'' J^ dt G{t)K^{at) . 
Some interesting examples from the above result are 



•''^ ^(x2-a2)(a;2+h2) LV ' V 7 I V V ; J 

r((p+i/+i)/2)r((p-i/+i)/2) p ftJ±j^+i tj.-u+i I -, . h^x 

r(77+Tl 2^*1 1^ 2 ' 2 '/''+ ^'^Z^J 

[Re M > -1, Im 6 = 0, Re (/i + 1) > jRe v\] 
[No. 2.1.9.21 Ref. [|, No. 2.16.21.1 Ref. @], 



(123) 



(124) 



(125) 



(126) 
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f dx 






^AF 



X t^Kf,{bx) - YJJITTJTj ^ 



'2-t^b^-Tiiy) B(^ + i,|-/,) ,F2(if^;l-z/,/i-| + l;4^ 
2-.-p-i„2.5M+.p(_^)5 (a* + i -/i - I) 1F2 (i±^; 1 + z., M + I + 1; 4^ 

2M-i„.-2M5-Mr (^ + 1) r (^ - I) 1F2 (i - m; 1 - M - 1, 1 + f - m; 4^) 



[Re ^ > -1/2, 6 > 0] 



[No. 3.16.1.6 Ref. [|, No. 2.16.3.7 Ref. g], 
and 



(127) 









-^'/^ A-^(bV^) = 2^"->-" ^ 



r(4^)r(^) 0F3 i 1 



i+i^ 1 I i^ — fJ- . a'b 
2 ' ' 2 ' 256 

2 y^v 2 yi'-^3i2! 2 



,2,1-4 ,y-M+i 



fcM- 



^\ 9 /^l 9 / 0^3 1 9 1 9 J 2 ' 256 



2_^_3.-i^2.5^+2.r(_^)r(-^ - ^i) ,f^ (i + 1^, i±^, i + ^; ^ 
23.-M-i5M-2.r(z.)r(z. - ^) 0F3 (1 - z., i±f^, 1 + 1^-9^ 

[|arg 6| < 7r/4] 

[No. 3.16.1.14 Ref. §, No. 2.16.8.9 Ref. §]. 
Specifically, if we put /i = in Eq. (127), then we find 



f dx 

J a 






KoibV^) = a-^Kl^, (^) 



or putting /i = 1, we get 



f dx 

Ja 



{xAr\/ x'^ —a^')^ -^{x— --J x^ —a^y^ —\ 
\/x^ —a^ 



Ki{by^) = a"-^ 



If we put ^ = and i/ = in (177), then we obtain 

r dx ^03 = 2[A:er2(5vW2) + kei\b./^)] 



(128) 



(129) 



(130) 



(131) 



We can develop the general result given by Eqs. (124) and (12J:) further by replacing v hy n + u and multiplying 
both equations by (— l)"/3"/n!. Then by summing from n = to oo and applying the multiplication theorem for the 
MacDonald function as given on p. 112 of Mangulis 0], we get 



f °° dx -^^ \{x + \/^^^^)''e-'3(-+VS^^^) + (2- _ Vx2 - a2).^e-;3(x-VJ^3^) 
2a'' J^dtG{t){l + 2(3/t)-'/^K^{a^/WTWt) ■ 



For the case of j/ = 0, Eq. (132) reduces to 



J^ dx ^gsL e-^^ cosh(/3Va;2 - a2) ^ J~ di G(i)i^o(av/^^T2^) 



(132) 



(133) 



In a similar manner, we can develop Eq. (121) further. This involves summing divergent integrals, the validity of 
which has been discussed in Sec. 3. Thus, by multiplying both sides of Eq. (121) by (— l)"(/3/2)2"+'^/n!r(n + i^ + I) 
and summing from n = to c», we eventually obtain after applying the multiplication theorem; 



J^ dx x'^/2(x - 2a)'^/2 J, ((3^x{x - 2a)) F{x) = ^^^^^ J^ dt e-"*G(i) x 

(^2 + ^2)-./2-l/4^^^^^^ (a^/WTW) , 



(134) 
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which is just a special case of Eq. JT^ ) after making the substitution y — x^l'^(x — 2a) ^/^. This, therefore, vahdates 
our technique of summing divergent integrals since Eq. ( [73[ ) was obtained using convergent integrals. 
To complete this appendix we consider integrals of the form 



/ = /q°° dx er]c{ax-^l'^)F{x) 
which can be transformed into 

An example arising from the above results is 



/q°° dx erfc(aa; ^1"^) {x^ + p 






23/2r(3/4)r(-i/) 



r ( — ij — l.\ P fl —i±l- 13 5. a'^ \ I p^"a^ p ( 1 _ 353. a"* 
^ V 4/ 2J^^3 \^4, 4 I 2' 4' 4' 64p2 ^ "r 4 2J^3 \^ 2 ' '^' 4' 4' 2' 64p2 



- p^-^/^g^ r(l/4-,.) 

2 6 r(i/4)r(-iy) 



f 3 1 

\^4' 4 



5 3 7. __a* 
4' 2 ' 4' 64p 



[Rei/ < -1/2, lmp = 0] 
[No. 2.1.5.1 Ref. [|, No. 2.12.9.4 Ref. §]. 



(135) 
(136) 



(137) 



VIII. APPENDIX C 



Here we establish the Feynman integral from its more general counterpart given as No. 2.2.6.1 in Prudnikov et 
al H and then develop further integral results from it. This integral, which appears so prominently in quantum 
electrodynamics, can be written generally |12| as 



Jo 



^a-l 



dx 



{l-x) 



0-1 



r(a)r(/3) 



{ax + b{l-x))"'+l^ T{a + P)a°'hP 



(138) 



where Re a > and Re /3 > 0. These conditions will apply throughout this appendix. To establish the above result, 
we shall show that 






(cx+d)-' 
2J 



(139) 



where Re a, /3 > and |arg(((i + cb)/{d + ca))\ < n. Eq. (139) is simply No. 2.2.6.1 from Ref. Q] with 7 replaced by 
—7. We shall primarily be interested in the case of Re 7 > even though the case of Re 7 < can be established by 
analytic continuation. 

With the aid of the integral representation of the gamma function Eq. (|139D can be written as 



r(7) 



(l-x) 



/3-1 



dtr- 



-ifx+g)t 



(140) 



where / = c{b — a) and g = ac+d. Expanding the part of the expo nential with fxt into a series and then interchanging 
the order of the integrations and the sum, we find that Eq. (140) yields 



(fc-Q)° 



Bic.,P) EZo 



(~l)"/" r(rt+a) r(a+/3) r(n+i) 
n! 



r(7) ^\^^I^J zL^n=0 n! r(a) r(n+a+/3) g' 

{b-a^+f^-^ Bia,f3) 2Fiia,r,a + (3;-f/g) , 



(141) 



which is just Eq. (139). To obtain the Feynman integral put a^O,b^l, c — a — b, d^b and 7 = a + /3. Then the 
2i^i hypergeometric function reduces to a iFq hypergeometric function, which in turn yields an algebraic function as 
given on p. 453 of Ref. Q. Thus one finally arrives at E q. (138). 

We shall consider the more general form of Eq. (138) with a + (3 replaced by 7 in the power of the denominator 
shortly. For now if we put a equal n + a with n a non-negative integer and multiply both sides by s"/n!, then we 
obtain after summing from n = to 00 and interchanging the order of the integration and summation 
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;o'^" (rJ+ 6(1-1))"+/^ e-/(-+^(^-)) = i3(a, /3) a-" 6-NFi (a; « + /?;./«) . 
On the other hand if we had multiphed both sides by (— l)"s"/n!r(n + 7 + 1), then wc would obtain 

-T/2-i(i^y3-i 



(142) 



/o dx 



J. 2 



{ax+b{l-x))''+f—i/^ '^■y I ^y ax+b{l-x) J " r(7+l) 

Alternatively, dropping the phase factor of (—1)" in the process yields 



X 



(143) 



J^dx 



-"-/2-1 (j^^-i 



U 2 



_ B{a£l 

{ax+b{l-x))'^+«-"'/^ -"T \^ ^V ax+b{l-x) ] r(7+l) 



w 



-iF2(a;7 + l,a + /3;f) 



(144) 



If we put 7 = a in Eq. 143, then by utilising tabulated results for 1F2 hypergeometric functions on p. 608 of Rcf. 
I we get 



r H-r ^ ('--X) J 9 / sx 

Jo "-^ (ax+b(l-x))''/^+f< "^" \ ^\ ax+b(l~x 



nl , x-'^-^ {l-x)^-' T I sx \ _ 2^+^-° .(^-°)/^r(/3) 



{ax+b{l-x))''/^+f< "" \."y ax+b{l~x) I ~ ai/2 bP 

{2aJa+p-l (2yf) SQ-/3-l,a+/3-2 {2^/^) - Jq+/3^2 (S^f) Sa-/3,a+/3-l (S^f)) , 

and for the case where a + (3 — 3/2, the above equation yields 



(145) 



/2-I (1^)1/2 



Jo "-^ (oa;+b(l-a;))(3°~i)/ 



Ja 2 



a I ^Y Q2:+6(l-2;) ) - h3/2-<, s=/2^ 



[r(2a - 1) cos(7ra) + S (2^, 2a - l)] 



(146) 



In Eq. (145), 5^,1^(2) represents the Lommel function while in Eq. (146), S{z,iy) represents the generalised Fresnel 
sine integral and Re a < 1/2. 



If we put 7 = a — 2 in Eq. (144), then we find using the results on p. 595 of Ref. g] that 
^1 , x-^- d-xf-' r . (0 I sx \ _ („_i),(g-^-i)/2r(ffl ^ 



h "-^ (aa;+b(l-x))'3+°/2+i ^a- 



ax-\-b{l—x) 



sl/2 



l,fsW+l)/2 



or if we put 7 = a — 1/2 and (3 ^ a, then we find 

1 ,./2-3/4(,_^).-l / / ,, \ _ ^^^ 

Jo ■^ (a£c+6(l-K))3°/2+i/4 ^a-1/2 I ^y ax+b{l-x) j ~ h° 50/2-1/4 

r(a + l/2)i?(a,a)/2_^/2(yT) . 
On the other hand, putting a = 1/2, 7 = 1 — c and (3 = c— 1/2 yields 



(147) 



,2c-l „-l/2 



(148) 



^'^/''-i (l-x)'' 



7rB(l/2,c-l/2) 



Jo '^^ (a2;+b(l-a:))3<=/2-i/2 -^l-c I ^Y'„2.+f,(l-2;) ) r(2~'c) sin(i^c)' ^ 



(149) 



the latter result being valid for < Re c < 2. 

We should also note that particular results for Eqs. ( |l43[ ) and (144) can be obtained by utilising the tabulated 
results of 1F2 hypergeometric functions in Prudnikov et al M with specific indices. For example, we find 



Jo ^^ iax+bil-x))^/f< '^1/4 [^y ax+b{l~-x) j " 6i/4s3/8 C (2,/^) 



(150) 
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Jo "^ ax+b{l-x} 



ax+b{l — x) 



s^(?Vl) 



(151) 



^^dx 



x-Uj-x)-^ 



{ax+b{l-x)y/ 



172 -^1 2 



ax+h{l — x) 



^(^MyD-AMVf)) 



and 



f^Wr ^'^^" /- , /^o / sx \ 

Jo "-^ (ax+b(l-x)y/i -'1/2 I ^Y aa;+6(l-£!;) I 



^/tt 6s'^ 



[cosh(2yf)-l] 



We add that many more results can be obtained by following t he a bove procedure. 

Now returning to the generalised Feynman integral, i.e. Eq. (138), we note that if we put a = n + a + 1, j3 
and multiply both sides by c^"+"/2^"+"ri!r(n + a + 1), then we obtain the following result after interchanj 
order of the sum from to oo and integration 



(152) 



(153) 



= n + /3 
^ing the 



(ax+b(l-x))l^+^ ^" \ ax+b(l~x) 1 ~ 2°a° + i bt> r(a+/3+l) 



Jo "-^ (a£!;+b(l- 

li'2\P, 2 '^- + -L'T6ZfcJ ' 

while if we multiply by (-l)"c2"+"/22"+"n!r(n + a + 1) instead, we get 

x"'^ (l-a,)g-°/^-i / c^x(l-x) \ _ c° r(/3) 



Jo ""^ (o2:+b(l-a:))3+i "^a I ax+b(l-x) j ~ 2°a° + i W* r(Q+/3+l) 

Some results emanating from the above are 



1 :,(T-3)/2 ^^_^^(^+i)/2 , cVa;(l-:r) \ __ 2-' r(7)r(7-l/2) 

Jo "-^ (a2;+f)(l-a;))-'+i -^7-3 I aa;+6(l-2;) I ~ a^/2-5/4 ;,t/2+3/4 r(27-2) 



-3/2 



c-y/a:;(l — a;) 
aj:+6(l— a:) 

-3/2 ViT^j ^ 8(7-l)\^6 -^T-l/a (iT^j 



(154) 



(155) 



(156) 



and 



Jo ^ {ax+bil-x))'-" ^1 I ax+b{l~x) 



^x(l~x) \ _ 2^ 
ac 



cosh 



\2Vab) 



J^dx 



{ax+b{l-x)y^ "'1 \ ax+b{l-x) 



-5/4 



^-^/^HV2 (^) 



In Eq. (158), H^(z) refers to a Struve function. In addition, by putting a — f3 — 1 in Eq. (|155|), one obtains 



Jo "^^ (a2;-|-f)(l-x))3+i -JP-i 



c7^(I3^ \ _ 2i-gc-V^V7 T f c \ 

ax+b{l-x) ) (afc)';V2+i/4 -'/3-1/2 1^2%/^^ 



(157) 



(158) 



(159) 



By utilising similar techniques to those described in this appendix, we can develop further results based on the 
Feynman integral such as 



Jo "-^ {ax+b{l-x))''+f^—i -'7 I QK+h(l-2;) ) ^ r(7+l)r(Q+/3)a° b'' ^ 



2^3 \^a, ^J 2 ' 2 ' '' 4ab 



(160) 



/o^ da; 2^ 



{ax+h{l~x)) 



F. aA 



r(a)r(l-a) 



6er 



( 



1. — z a;(l — a;) 
'^' 2' 16(aa;+b(l-a:))^ 



Re a < 1 



(161) 
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A similar result to Eq. ( 161 ) can be found for the Kelvin function, bei{z), by using No. 8.564 from Gradshteyn and 
Ryzhik [|l|. Finally, Eq. ( 160 ) yields interesting results of its own, two of which are 



and 



a-" b-P iFi (a; a + /3; y^) i^i (a; a + /?; -y^) , 



(162) 



/o da: 



{ax+b(l-x)Y/^ 



ha 



2y/cxil-x) \ _ 2'i''-^c^/^ 
ax-{-b{l — x) 



B(a,a + l/2)/„_i 



/4 



ai/2 6 



iZ&J ^a-3/4 



B(a + l/4,a + 3/4) x 



(163) 
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